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Abstract 

We confront the hypothesis of chemical freeze-out in ultra-relativistic heavy-ion 
collisions with the hypothesis of large modifications of hadron masses in nuclear 
medium. We find that the thermal-model predictions for the ratios of particle mul- 
tiplicities are sensitive to the values of in- medium hadronic masses. In particular, 
Q\ , the 7r + jp ratio decreases by 35% when the masses of all hadrons (except for pseudo- 

Goldstone bosons) are scaled down by 30%. 

PACS: 25.75.Dw, 21.65.+f, 14.40.-n 

Recent theoretical studies [1,2] indicate that hadron yields and ratios in ultra- 
relativistic heavy-ion collisions agree with predictions of a simple thermal 
model. According to this model all measured abundances of hadrons are con- 
sistent with the assumption of complete thermalization of hadronic matter at 
a temperature T chem , baryon chemical potential /i^ em , strangeness chemical 
potential ^ hem1 and isospin chemical potential /i^ em . Thermodynamic param- 
eters obtained from this type of analysis characterize the point in the evolution 
of the system when "primordial" hadron content is established. One refers to 
this point as to the chemical freeze-out. At this stage the system is a mix- 
ture of stable particles (pions, kaons, nucleons,...) and resonances (p, u, A, ...). 
In the subsequent evolution the resonances decay, contributing to the final 
(observed) multiplicities of stable particles. 

It has been shown that the chemical freeze-out parameters at CERN/SPS, 
BNL/AGS and GSI/SIS energies all correspond to a unique value of the en- 
ergy per hadron [3]. Moreover, statistical models [4] are able to reproduce 
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the particle yields in e + e~ collisions. A striking observation is that at very 
high energies the temperature T c h em ~ 170MeV turns out to be the same for 
both elementary and nuclear collisions, although the final-state hadronic in- 
teractions are completely absent in the former case. This may indicate that 
chemical equilibrium is pre- established by the hadronization process [4-6]. 

In recent studies one distinguishes the chemical freeze-out from the thermal 
or kinetic freeze-out [5]. The latter is defined as the decoupling of strongly 
interacting matter produced in high-energy nuclear collisions into a system of 
essentially free-streaming particles. After the thermal freeze-out the hadrons 
practically stop to interact and travel freely to detectors. Chemical freeze-out 
is expected to occur at the same time or before the thermal freeze-out [5]. For 
Pb+Pb collisions at CERN/SPS energies one finds that the chemical freeze- 
out point occurs significantly earlier than the thermal point. This is indicated 
by the measurements of the hadron momentum spectra as well as two-particle 
momentum correlations [7], which show that the thermal freeze-out temper- 
ature T^erm ~ 100 MeV is substantially lower than T c ^ em ~ 170 MeV. In 
addition, T c h em is very close to the expected critical value for the deconfine- 
ment/hadronization phase transition, therefore the chemical composition of a 
hadronic system, according to the presented scenario, should be established 
shortly after hadronization. 

Since the chemical freeze-out occurs at an early stage of the evolution of 
the system, where temperatures and densities are very high, we expect that 
hadronic properties, such as masses, widths, or coupling constants are in- 
fluenced strongly by the medium. Such modifications are predicted in many 
theoretical calculations [8-11]. Moreover, they seem much desired [12,13] in 
view of the low-mass dilepton enhancement observed in the CERES [14] and 
HELIOS [15] experiments. 

In order to study how the chemical freeze-out parameters are affected by the 
in-medium change of the hadron masses, we calculate the particle densities 
from the standard ideal-gas equilibrium expression 
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where gi is the spin degeneracy factor of the ith hadron, Bi, Si, Ii are its 
baryon number, strangeness, and the third component of isospin, and E* = 



p 2 + (m*) is its energy. The latter explicitly depends on the in-medium 
hadron mass m*. 

In the thermal-model fits [1,2] one uses Eq. (1) with vacuum values of hadron 
masses m* = m^, and, in addition, applies finite-size and excluded- volume 



corrections. They account for the finite size of the nuclear system and the 
finite volume occupied by the individual hadrons. The main effect of such 
corrections is a reduction of the absolute yields of particles with the particle 
ratios remaining close to predictions of the ideal-gas approach [2]. For Pb+Pb 
collisions at CERN/SPS energies the predictions of the thermal model are [2]: 

T chem = 168 MeV, ^ hem = 266MeV, 

(2) 
»Lm = 71MeV, /4 em = -5MeV. 



In our study we use these values in Eq. (1) and calculate the particle densities 
Hi. We take into account all meson resonances with masses smaller than 1.28 
GeV and baryon resonances with masses smaller than 1.45 GeV. Decays of 
resonances contribute to the final (measured) yields of stable hadrons. This 
is an important effect [1]. The needed branching ratios for the decays of reso- 
nances are taken from experiment [16]. In this paper, for clarity, we neglect the 
finite-size and the excluded-volume corrections. As mentioned above, they do 
not change significantly the ratios of particles. Indeed, within our simplified 
approach we reproduce the numbers of Ref. [2] at the level of 15%. 

In principle, the in-medium masses of all hadrons may behave differently. 
For practical reasons we perform our calculation with the meson and baryon 
masses rescaled by the two universal parameters, xm and xb, namely 

m* M = x M m M , m* B =x B m B . (3) 



In this way we change the masses of all hadrons except for pseudo-Goldstone 
bosons, i.e., pions, kaons and the eta, whose masses are kept fixed at vacuum 
values. In fact, results of many model calculations show stability of the pion 
mass against the change of temperature and density up to the point where 
the chiral phase transition occurs [17]. The branching ratios are kept at the 
vacuum values. 

In Fig. 1 we show the 7r + /p ratio calculated for different values of the parame- 
ters Xm and xb- The solid line represents the case when the meson and baryon 
masses are rescaled the same way, xm = %b = x (BR-scaling [8]), the dashed 
line corresponds to the case when only the baryon masses are changed, xm = 1 
and xb = x, and the dotted line shows the case xm = x and x B — 1- In all 
three cases the values of xm and x B are restricted to the reasonable range of 
0.6 <x< 1.1 

Our results presented in Fig. 1 show that large mass modifications lead to a 
substantial change of the 7r + /p ratio. The reduction of masses by 20% changes 
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Fig. 1. The ir + /p ratio plotted as a function of the scaling parameter x. Solid line: 
all hadron masses (except for Goldstone bosons) are scaled with xm = xb = x. 
Dashed line: only baryon masses are scaled xu = 1> %B = %■ Dotted line: only 
meson masses (except for Goldstone bosons) are scaled xm = x, xb = 1- scaled 

the 7r + /p ratio at least by 30%, which occurs in the case where xm = %b- In 
the case xm = 1 a decrease of baryon masses by 20% reduces the 7r + /p ratio 
by 50%. On the other hand, in the case of constant baryon masses a decrease 
of meson masses by 20% enhances the 7i + /p ratio by 60%. The displayed 
behavior can be easily understood in qualitative terms. The strong increase of 
the dotted curve in Fig. 1, as the meson masses are decreased, is mostly caused 
by a larger population of the rho and omega mesons at the chemical freeze- 
out point. Subsequent decays of these mesons produce more pions, raising the 
7i + /p ratio. We note that the rho and omega bring about half of the effect 
shown in Fig. 1. The other half comes from heavier resonances. Similarly, for 
the dashed curve, a lower mass of baryons results in the larger population 
of protons at the chemical freeze-out point, thus lowering the 7i + /p ratio. A 
universal scaling of the meson and baryon masses (solid line) partially cancels 
the above-described effects. Still, a significant effect remains. Other particle 
ratios are less sensitive to the mass modifications. 

To conclude, we note that the thermal fit analysis leaves little room for mod- 
ifications of hadron masses larger than, say, 10-20%. The quite impressive 
agreement with data reached in [2] would be deteriorated with significant 
mass modifications. However, it seems premature to jump to a general con- 
clusion that masses of hadrons cannot significantly change in hot and dense 



medium. One should bare in mind that at present we do not understand in 
sufficient detail the mechanisms of hadron production and evolution of the 
system created in heavy-ion collisions, and the appealing simplicity and nu- 
merical success of the thermal model may be misleading. One cannot exclude 
the possibility that a more elaborate or altogether different treatment will 
leave room for a substantial modification of hadron masses. In order to shed 
more light on this issue it would be useful to incorporate scaling of hadron 
masses, as well as other hadronic parameters, in various existing approaches 
to heavy-ion collisions. 
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